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The 5 = 1/2 Ising chain with uniform nearest-neighbor and next-nearest-neighbor coupling is used to construct 
a system of floating particles characterized by motifs of up to six consecutive local spins. The spin couplings 
cause the assembly of particles which, in turn, remain free of interaction energies even at high density. All mi- 
crostates are configurations of particles from one of three different sets, excited from pseudo-vacua associated 
with ground states of periodicities one, two, and four. The motifs of particles and elements of pseudo-vacuum 
interlink in two shared site variables. The statistical interaction between particles is encoded in a generalized 
Pauli principle, describing how the placement of one particle modifies the options for placing further particles. 
In the statistical mechanical analysis arbitrary energies can be assigned to all particle species. The entropy is 
a function of the particle populations. The statistical interaction specifications are transparently built into that 
expression. The energies and structures of the particles alone govern the ordering at low temperature. Under 
special circumstances the particles can be replaced by more fundamental particles with shorter motifs that in- 
terlink in only one shared site variable. Structures emerge from interactions on two levels: particles with shapes 
from coupled spins and long-range ordering tendencies from statistically interacting particles with shapes. 
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1. Introduction 

Condensed matter is an aggregate of interacting particles. The interactions operate in hierarchies 
regarding strength and range. Atomic nuclei are composed of strongly coupled protons and neutrons. 
Electrons are bound to atomic nuclei by electromagnetic interactions of widely varying strength. Inner 
electrons are tightly bound to the nucleus and form ion cores. Some outer electrons may be traded or 
shared between neighboring ion cores in ionic or valence bonds, respectively. In metallic bonds some 
outer electrons are mobilized. 

Sorting out the diverse, complex, and interrelated phenomena is challenging. The common strategy of 
many approaches is to transform specific aspects of the strongly interacting ion cores and electrons into 
more weakly interacting collective modes. The latter share many attributes with fundamental particles 
including energy-momentum relations, spin, and exclusion statistics. They scatter off each other elasti- 
cally or inelastically, form bound states, or decay into other modes. The goal of transforming the strongly 
coupled constituent particles of condensed matter into collective modes that behave like free particles is 
elusive except under idealized circumstances related to dynamic or kinematic restrictions. 

Harmonic lattice vibrations have linear equations of motion. The collective modes exist in superpo- 
sitions of infinite lifetimes without scattering. Phonons have bosonic statistics. Linear combinations of 
atomic orbitals produce fermionic counterparts: band electrons |lj]. In a model that constrains the kine- 
matics of collective modes to one dimension and limits their dynamics to elastic two-body scattering, 
the momenta and energies are conserved individually. The collective modes are particles with exclu- 
sion statistics determined by the (factorizing) S-matrix, analyzed via Bethe ansatz Lattice degrees 
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of freedom that are coupled by commuting operators are static. They can be assembled into particles 
that are free and floating. These particles have definite energies independent of the neighborhood. Their 
mutual exclusion statistics are exotic 

This project imposes the kinematic constraint of one dimension and the dynamic constraint of com- 
muting operators. It aims at shedding light on the assembly of structures from interactions and on the 
emergence of order from these structures within that limited realm. We expect that the results will open 
doors for situations where the kinematic or dynamic constraints are relaxed. The Ising chain for commut- 



ing spin operators Sf 



with nearest-neighbor (nn) coupling /, next-nearest-neighbor (nnn) coupling 



L, magnetic field h, and periodic boundary conditions \ 9}, 



^ r 



(1) 



is a suitable starting point from which it is possible to develop the methodology. The spectrum consists 
of product states, \a\---aN) wither; The notation |(Ti ■■■(Tjv)p refers to sets of p product states (of 
periodicity p) that transform into each other via translations. 

The zero-temperature phase diagram aXh-Q features three phases with periodicities p- 1,2,4 (see 
figurejl}- This includes phase Oi with nn spins aligned, phase $2 with nn spins anti-aligned, and phase <1>4 
(at L > \ \]\) with nn spins alternatingly aligned and anti-aligned. The nnn spins are uniformly aligned in 
phases (Pi, <1>2 and uniformly anti-aligned in phase $4. At ^ and / > 0, phases <1>2 and <1>4 persist, phase 
$1 is split up into two phases <l>i+, and two new plateau phases $3+ with periodicity p-i are stabilized^ 
The hi^Q phase diagram at L = features phases <1>2 in the sector J >\h\ and <l>i+ at % in the other 
sectors. 
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Figure 1.7 = phase diagram (a) for positive nn coupling / and nnn coupling L of either sign in the 
parameter plane of scaled nnn-coupling LI J and scaled magnetic field hi J, (b) for h = Q in the {J,L)- 
plane, and (c) for L = in the (/, h)-plane. The text refers to phases and regions in parameter space by the 
same names. 



In section[2]we configure the physical vacua O2, and <1>4 as the pseudo-vacua of sets of statistically 
interacting particles. Salient features of the statistical mechanics of these particles are highlighted in 
section [3] (entropy, ordering tendencies) and section |4] (populations in competition). The emergence of 
structures from interactions on two levels is further discussed in section [s] 



-'The hi^Q phase diagram is much simpler at /< 0. Only three of the phases exist: <1)4 in the sector L/|/| > + 1/2 and<I)i+ at 
?! S in the remaining sector of the (L/|/|, /!/|/|)-plane. 
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2. Motifs, species, categories 

Motif is a term borrowed from music, literature, and visual arts, where it refers to fragments, themes, 
or patterns. Now well established in the natural sciences, it is used in biochemistry, for example, to de- 
scribe patterns of nucleotides (codons) in DNA sequences. Motifs as a representation of particles in many- 
body systems are common in statistical mechanics . In the present context, motifs are patterns of 
consecutive site variables ct; in Ising product states. Each motif characterizes either a particle of a partic- 
ular species or an element of pseudo-vacuum. Our goal is to find motifs representing particles that have 
a definite energy irrespective of their location on the lattice relative to other particles. This requirement 
limits the number of ways particles can be identified. Every product eigenstate |fJi • • • ctjv) thus becomes 
a string of motifs. 

Successive motifs interUnk in a characteristic manner that depends on the range of interaction, (i) 
The h-tevm in Hamiltonian (1) has zero range. It permits the use of motifs that interlink by sharing no 
site variables (e.g. | -i- 1 = ff )■ All on-site energies are additive, (ii) The /-term has range one. It demands the 
use of motifs that interlink in one shared site variable (e.g. H + f I = I f I ). This makes all n n-bond energies 
additive and the additivity of on-site energies can be maintained by specific rules, (iii) The L-term has 
range two, which requires that motifs interlink in two shared site variables (e.g. nf + ff| = nf|). This 
guarantees that all nnn-bond energies are additive without jeopardizing the additivity of nn-bond and 
on-site energies. Longer-range couplings demand a more extensive overlap of motifs and more elaborate 
rules for ensuring additivity of energy contributions from all Hamiltonian terms. 

The search for motifs of a set of free particles that are excited from a given pseudo-vacuum is guided 
by the additional optimization criteria that aim for the fewest and shortest motifs. In the context of a spin- 
1 Ising chain with nn-coupling we identified sets of six particles excited from twofold pseudo-vacua and 
sets of seven particles from non-degenerate pseudo-vacua |8]. Here we use the same strategy to identify 
three sets of particles that generate the full spectrum of (1) from pseudo-vacua | H •••)!. IIUI •••)2. 
I 1 1 1 U I U • • • >4, associated with the ground states at h-0. 

The taxonomy of particles defined by their motifs involves structures as sorted into species and func- 
tions as emerging from categories. These features were the focus of reference Q]. In the present context 
the species will be very different but the categories will remain the same. We shall again encounter 
compacts, hosts, tags, hybrids, and no further categories. Compacts and hosts float in segments of pseudo- 
vacuum, tags are located inside hosts, and hybrids are tags with hosting capability. Particles from the 
same category but with different structures may collectively allow the emergence of new functions. 

The number of product eigenstates that contain specific numbers {N,n] of particles from all species 
of a given set is expressible by a multiplicity function W{{Nm]]- Its general structure, developed in the 
context of reference |[8(], remains operational without modification: 

J" J" m=l ^ ' m=l 

M 

dm — Afn~ ^ Smm'i^m' ~ ^ mm')> (2b) 
m'=l 

where ripy is the multiplicity of the pseudo-vacuum, the Am are capacity constants, the are size con- 
stants, and the gmm' are statistical interaction coefficients. The generalized Pauli principle proposed by 
Haldane ] 11] is encoded in (2b) with dm counting the number of open slots for particles of species m in the 
presence of N^' particles from any species m', thus encapsulating the essence of statistical interaction. 
All product states with particle content [Nm] have energy 

M 

E[{Nm]]^Ep^+ NmCm, (3) 
m-1 

where Cm is the energy of particles from species m relative to the pseudo-vacuum, which has (absolute) 
energy £pv. 
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2.1. Particles generated from | titi •••>2 

The physical vacuum in region <I>2 is the twofold Neel state | fif ■■■)2, here selected as the pseudo- 
vacuum for M - 4 species of particles with specifications as compiled in table [l] . The permissible con- 
figurations of particles from these species generate the complete spectrum of The motifs interlink as 
illustrated in figure [2] Particles m - 3,4 are hosts and particles m - 1,2 tags. Hosts can be placed into 
segments of pseudo-vacuum and tags inside hosts. In this instance, hosts 3 accommodate tags 1 only and 
hosts 4 tags 2 only. 



Table 1. Specifications ofM = 4species of particles excited from the Neel state (npv = 2) | ][] •••>2: motif, 
category, species, energy (relative to pseudo-vacuum), spin, capacity constants, size constants (left), and 
statistical interaction coefficients (right). Segments of ( vacuum elements, t i t , It 1. have energy l?[L-J)li. 
AthtO the entries of e„i must be amended by -s^h. 



motif 


category 
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Cm 
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host 
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host 


4 
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N-1 
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8mm' 
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1 








-1 
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-1 
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1 
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1 


3 
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2 


2 


2 



The energies of hosts and tags are calculated differently. In hosts (or elements of vacuum) we count 
interior nn bonds fuUy but the outermost nn bonds only half. In tags we do not count the nn bond on the 
left but do count the one on the right fully. All nnn bonds are counted fully in each particle (or element 
of vacuum). These rules are implemented in the entries for In analogous manner a spin Sm can be 
assigned to particles from each species. 

Tags have vanishing Am ill. Open slots for tags are created by the placement of hosts. The number 
of tags that can latch on to a given host is only limited by the space to which the latter can expand in the 
process. In one product state tags (of one species) are lined up around the chain with no need of a host. 

The gmm' from table[l]as used in l [2b) may be interpreted as follows: (i) Adding a tag does not change 
the capacity for further tags. Any slot taken by a tag opens up exactly one slot for a tag of the same kind, 
(ii) In the process of adding a tag its host must expand, which reduces the capacity for further hosts of 
either kind in the same manner, (iii) Adding a host increases the capacity for tags that it can host and has 
no effect on the capacity for tags it cannot host, (iv) Adding a host of any kind diminishes the capacity 
for further hosts of the same kind more strongly than for hosts of the other kind. Different hosts can be 
interlinked directly, whereas identical hosts must be spaced by at least one element of vacuum. 

The two vectors of phase <1>2 are the pseudo-vacuum of the particles considered here. The two vectors 
of phase Oi are each a solid of N close-packed tags of one kind. By contrast, the four vectors of phase <1>4 
are each a solid of ^ A? hosts in an alternating pattern. 
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Figure 2. Interlinking elements of pseudo-vacuum and particles from table [T] Some nn bonds (u) are 
counted half, others (v) full. All nnn bonds (a) are counted fuUy. 
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For the case of vanishing nnn-coupling (L - 0) the four species of particles from table[T]can be replaced 
by the two species in table|2] The motifs are shorter and interlink in one shared site variable. All nn bonds 
are now counted fuUyand aU sites half in each motif. The new particles (named solitons) are well known 
from previous work J^Ji3\. Their statistical interaction is semionic, an attribute shared with the spinous 
identified in the XX Z and Haldane-Shastry models |13, 14]. In the taxonomy of reference |8] they are 
compacts. In figure[l](c) phase ^2 is the pseudo-vacuum of solitons whereas phase O1+ (<l>i_) is a solid of 
N spin-up (spin-down) solitons. 



Table 2. Specifications of M = 2 species of particles excited from the Neel state for L = (n 



pv 



2) 



I till •••>2- Segments of £ vacuum elements, tl.lt, liave energy -£Jli. At htO the entries of must 
be amended by -Smh. 



motif cat. m 




1 1 comp. -1- 
i i comp. - 


2^22 ^ 
J 1 iV-l , 
222^ 



Emm' 


-1- 


-1- 


1 1 

2 2 

1 1 

2 2 



The motifs of solitons are not fragments of the original motifs. They interlink differently. At L it we 
have two spin-up particles (tag 1 and host 3) and two spin-down particles (tag 2 and host 4). At L = 
the energies of both spin-up (or spin-down) particles become equal. In this case we can get away more 
economically with one spin-up and one spin-down soliton. Returning to L ^t 0, the energy of a soliton 
depends on its position relative to other solitons. The interaction energy is eliminated by switching back 
to the extended set of four species. 



2.2. Particles generated from 1 1 1 - - - ) 1 

Phase <l>i comprises the twofold spin-polarized state | \ \\ ■■■)\, \ \\ \ ■■■)\. Here we adopt the first 
vector as the pseudo-vacuum for M = 5 species of particles with specifications compiled in table[3] Three 
categories are represented. Host 2 accommodates only tag 4 whereas host 1 accommodates both tags 
and the hybrid. The hybrid, in turn, is capable of hosting both tags. The energies of hosts and tags are 



Table 3. Specifications of M = 5 species of particles excited from the spin-polarized state («pv - 1) 
I n ■ • ■ > 1 • Segments of ( vacuum elements, ]]], have energy £{L + J)I4:. AthtO the entries of must be 
amended by -Smh. 
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m 


Cm 








imn 
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2 
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3 
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calculated as in section [2A] The rules for hybrids are the same as those for tags. The pseudo-vacuum is 
spin-polarized and Sm is not a spin in the usual sense. It enables us to write the magnetic-field contribu- 
tion to Cm in the form -Smh. 

Various combinations of hosts, tags, and hybrids are illustrated in figure[3] In the search for particles 
that are free of interaction energies, attention had to be paid to the requirement that the implantation of 
a tag or hybrid into a host (or a tag into hybrid) leaves the sums of aligned and anti-aligned nnn bonds 
invariant. 



fiti = tAtitt 
itit = ttititt 

fiti = ufitu 



Figure 3. Interlinking particles from table [3] The number of nnn bonds is preserved and the number of 
nn bonds is reduced by one. The shared bond is marked by a or v. 

The hosting capabilities of particles m' - 1,2,5 are encoded in negative statistical interaction coeffi- 
cients gmm'- In all instances except one, we have gmm' - The fact that host 1 has two interior slots to 
accommodate tag 4 requires that g4i - -2. Tag 3 has zero energy, a consequence of our choice of pseudo- 
vacuum. The physical vacuum in region <l>i comprises the pseudo-vacuum | ff • • • ) i and a solid of tags 3, 
I ii ■ ■ ■ ) 1 . In region <I>2 the physical vacuum is a solid of negative-energy tags 4 (two vectors) and in region 
<I>4 a soUd of negative-energy hosts 1 (four vectors). 

For the case L = we can again get away with fewer particles that have shorter motifs and interlink 
with shorter overlap: one host and one tag as listed in table H) In figure [l](c) phase $1+ is the pseudo- 
vacuum, phase Oi- is a sohd of tags, and phase $2 is a soUd of hosts. The two species of table|4]are free of 
nn interaction energies. The nnn interaction energies between them can be eliminated if we allow them 
to assemble into the five species of table[3] 



1-1-2 ttUtt + 

1-1-3 ttAtt + 

1+4 tAitt + 

1+5 ttAtt + 



ttitt = ttUttitt 

Ai = ttuitt 

■ At! = tAtutt 

. i;it = tti^titt 

Atu = ttAtiitt 



2-1-4 tt\tt + 
5-1-3 utA + 
5+4 u/u + 



Table 4. Specifications of M = 2 species of particles excited from the spin-polarized state for L = («pv = 1) 
I It •••)!. Segments of £ vacuum elements, ]], have energy £JI4. Mh + Q the entries of must be 
amended by -Smh. 



motif cat. m 




in host H 
U tag T 


-J -\ N-l I 
-I 1 



Smm' 


H T 


H 
T 


2 1 
-1 



2.3. Particles generated from 1 1 1 i i 1 1 • • • >4 

In region (S>4 the fourfold state | If I U I • • • )4 is the physical vacuum. Our search for free particles that 
generate the entire spectrum from this state configured as pseudo-vacuum has produced M- 4 compacts 
with motifs and specifications compiled in table[5] AU particles or elements of pseudo-vacuum comprise 
two nn bonds and one nnn bond. They again interlink by sharing one nn bond. Their energy content 
consists of one half of each nn bond energy plus the fuU nnn bond energy. 

The assignment of a spin ± ^ to the four particles is based on the following reasoning. Take two inter- 
linked elements of pseudo-vacuum with zero spin, ff U or I U I. add one of the four particles, and check 
the spin of the resulting entity. For example, H U + 1 1 1 - 1 1 1 U and H U + I U - A U I produce entities 
with spin -1- j . 
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Table 5. Specifications of M = 4 species of particles excited from the state («pv = 4) Ittlltt •••>4- Seg- 
ments of ^ vacuum elements ttiitlLUtiltt have energy -£LI4. At ^ the entries of must be 
amended by -Sm h. 



motif category m 




1 1 1 compact 1 
i i i compact 2 
1 i 1 compact 3 
i 1 i compact 4 


i(2L + /) +i ^ 1 
i(2L + /) -i ^ 1 
i(2L-/) +i iv±i 1 
\[2L-J) -\ ^ 1 
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1 
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4 
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1 


1 


3 




4 


4 


4 


4 



Four coefficients gmm' are negative. Here we do not have the standard host-tag scenario described 
previously. To understand the variant scenarios we note that each motif (four elements of vacuum and 
four particle species) can be followed only by two out of eight motifs. Only the motifs of particles 1,2 
can follow themselves. Thus adding a particle m' = 1 or m' = 2 merely reduces the total number of open 
slots for further particles of any species. This explains that all coefficients in the first two columns are 
equal and positive. The coefficients in columns m' = 3,4 reflect stronger exclusion in some instances and 
accommodation in others, involving two distinct mechanisms. 

(i) We note that particles 3,4 accommodate each other mutually. The addition of a particle m' - 3 
opens up a slot for a particle m- 4 and vice versa. This mutual accommodation is reflected in the 
coefficients g34 - g4s - -1/4. At the same time the addition of a particle m' = 3 or m' = 4 closes 
down a slot for further particles of the same species. This accounts for the more strongly positive 
coefficients gas - g44 = 3/4. 

(ii) The sequence between any motif and the motif of particle m-l can always be shortened by the in- 
sertion of a particle m' -3 wherever it fits. By contrast, the insertion of a particle m' - A wherever 
it fits will lengthen that sequence. The presence of a particle m' -3 thus increases the capacity of 
the system for particles m-l and the presence of a particle m' - 4 has the opposite effect. This 
is reflected in the coefficients gis - -1/4, gi4 - 3/4. Analogous reasoning explains the tabulated 
values of g24 and g23. 

The ferromagnetic (FM) phase <l>i consists of two states with broken spin-flip symmetry, one being a 
solid of particles 1, the other a solid of particles 2. The antiferromagnetic (AFM) phase <1>2 also consists 
of two states, each solid composed of particles 3,4 in an alternating sequence with broken translational 
symmetry. FM particles 1,2 and AFM particles 3,4 both have spin Sm = ± | . By interlinking differently, the 
former produce a uniform magnetization in phase Oi and the latter a staggered magnetization in phase 
O2. 



3. Entropy landscapes 

The statistically interacting particles from the three sets identified in section[2]have definite shapes 
and energies. These floating objects are assembled from localized spins by the nn and nnn couplings of 
Hamiltonian Q. The Ising chain is an open system of particles with energies depending on J,L,h. 
Here we abandon the Ising context and focus entirely on particles with interlinking motifs. 

The statistical interaction between the particles depends on their shapes and on the nature of the 
pseudo-vacuum, but not on the particle energies Cm- It is instructive to explore the effects of statistical 
interactions produced by particular shapes in a setting where particle energies do not factor in. To this 
end we consider the configurational entropy as derived from the multiplicity expression for N, N,„ » 
IviaS^fcelnW: 

M 

s[{N„,}) = fcfi E K^'" + 1" (^"^ + ^'») - i"^'" - ' (4^^ 

m=l ^ ' 
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M 

Ym^Am- ^ gmm'Nm'- (4b) 
m'=l 

The functional dependence of S on the populations Nm of particle species from a given set determines an 
entropy landscape shaped by the statistical interactions alone. 

If only one species is present, the statistical interaction reduces to an exclusion principle. The function 
S{Nm) vanishes identically for any species of tags or hybrids, which have Am - and gmm - 0. These 
particles can only exist inside hosts. Compacts and hosts exist in segments of pseudo-vacuum. They have 
Am oc N and gmm > 0. The entropy S(Nm] of a single species of hosts or compacts is nonzero for < 
Nm < Am/ gmm and zero at the endpoints. 

The signature of the statistical interaction between any two species m, m' is best visualized in a re- 
duced entropy landscape, in contour plots of the entropy per site, S{Nm,Nfy,i) - S{NmlN,Nm' IN) IN with 
the populations of all other species suppressed. All thermodynamic processes described in the following 
are understood to be quasi-static and to be implemented in an open system^ Equilibration would be 
problematic in closed, one-dimensional systems of particles from more than one species. 

3.1. Hosts, tags, hybrids 

Beginning with the four species from section lZTI we consider the entropy landscapes pertaining to the 
two hosts, a host and the tag it does or does not accommodate, and the two tags. In each instance we vary 
the population densities of two species over the permissible range while keeping the population densities 
of the other two species constant at very low values. 

The entropy landscape of the two hosts as shown in figure |4] (a) has borders of quadrilateral shape. 
The four corners correspond to pure phases with vanishing entropy. Hosts interlink with elements of 
pseudo-vacuum in multiple configurations. This explains the nonzero entropy along the two sides on the 
axes. The other two sides represent states that are crowded with hosts of both species. If - N4 the 
state of maximum population density is unique. Hosts 3,4 interlink directly. Hosts 3,3 or 4,4 are spaced 
by at least one element of pseudo-vacuum. The corner states '1>3+ contain close-packed hosts from a 
single species. The state <1>4 contains jN close-packed hosts from both species arrayed in an alternating 
sequence. Heading from point <1>3+ toward point <1>4 involves the repeated replacement of one host 3 by 
three hosts 4. 

A qualitatively different entropy landscape pertains to host 3 and tag 1 (hosted by 3) as shown in 
figure|4](b). The borders are of triangular shape. Hosts 3 alone (vertical leg) generate entropy but tags 1 
alone (horizontal leg) do not. All tags 1 are arrayed uniformly inside hosts 3. Only the expanded hosts of 
multiple sizes have positional disorder, not the tags inside. As the hosts disappear, so does the entropy. 
The tags contribute to the entropy only indirectly by expanding hosts to different sizes. Close-packed 
configurations of hosts stuffed with tags are represented by points on the hjrpotenuse. Near the middle 
the entropy is largest, generated by close-packed hosts expanded to many different sizes. Near one end, 
we have many more hosts but with only few small sizes represented. Near the other end, we have few 
hosts expanded to large sizes. Either trend reduces the entropy. 

The entropy landscape is yet different for host 4 and tag 1 (not hosted by 4) as shown in figure|4](c). 
AU tags 1 are confined to very few hosts 3, incapable of producing any significant entropy by themselves 
(horizontal leg). Almost all hosts 4 have the minimal size due to the near absence of tags 2. These hosts 
4 produce positional entropy along the vertical leg. The hypotenuse describes configurations of near 
uniform segments of hosts 4 separated by rare hosts 3 filled to various sizes with uniform arrays of tags 1. 
This arrangement has very low entropy. Starting with any density of hosts 4, the entropy always decreases 
when we add tags 1. AU tags 1 added replace elements of pseudo-vacuum and thus reduce the options for 
positioning hosts 4. The tags themselves are lumped together inside few hosts 3. 

The two tags 1 and 2 have an entropy landscape (not shown) that is also triangular but very flat near 
zero entropy. The two species of tags do not mix. They are arrayed uniformly in separate hosts, present 
only in very low densities. 

The entropy landscapes of the five species of hosts, tags, and hybrids from section [2^ are all of trian- 
gular shape. The hybrid 5 gives rise to one new feature, e.g. in combination with tag 3 as illustrated in 

^The Ising model Q is one particular realization with specific particle energies £,„ ■ 
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Figure 4. (Color online) Entropy per site, Slk^, versus population densities, Nm , , of two species from 
(a)-(c) table[l]and (d) from table[3] The population densities of the other species are kept at 10~^. The 
contours are at (a) 0.071f, (b) 0.068f, (c) 0.039^, (d) 0.058^, £ = 1 6. 



figure |4]((i). Both species exist inside hosts 1 or 2. Tags 3 alone or hybrids 5 alone produce no significant 
entropy for reasons already stated. However, the two species can coexist inside the same host in many 
different configurations. Almost all entropy is now generated inside the few hosts available. 

In surveying the four panels of figure |4] we see that the entropy maximum is realized under diverse 
circumstances, with both species or only one species present at intermediate density, or with two species 
in close-packed configurations. As expected, maximum positional ordering of the particles (zero entropy) 
occurs if no particles are present or if the system is close-packed with one species. Close-packing with two 
species is also found to produce maximum positional order either through alternate stacking [panel (a)] 
or through segregation [panel (c)]. 

3.2. Compacts 

All four species of compacts from section lOl interlink with elements of pseudo-vacuum but only se- 
lectively with themselves or each other. This produces two features in their entropy landscapes not seen 
in the previous cases. Compacts 1 and 2, which interlink with themselves but not with each other or 
with any other compact, generate a triangular entropy landscape as shown in figure[S](a). The entropy is 
nonzero along each leg, where compacts from one species mbc with elements of pseudo-vacuum. The en- 
tropy is zero along the hypotenuse, where compacts are close-packed. The system has the highest capacity 
for compacts 1 and 2 if they are segregated. 

The most remarkable entropy landscape pertains to compacts 3 and 4 as shown in figure[5](b). These 
two compacts interlink with each other but not with themselves or with any other compact. The border is 
a quadrilateral as already seen in figure|4](a) but with the sides away from the axes now slanted positively. 
If we start fUling the system with compacts 3, the entropy first increases from zero and then decreases 
again down to zero when capacity is reached at N3- |w in a uniformly stacked array. Interestingly, we 
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Figure 5. (Color online) Entropy per site, SI k^, versus population densities, Nm , N^i , of two species from 
table[5] The population densities of the other species are kept at 10~^. The contours are at (a) 0.070f , (b) 
0.071^, ^= 1,..., 6. 



do not have to remove any compacts 3 to make space for compacts 4. In fact, for every three compacts 4 
added we can add one more compact 3. The entropy rises in the process and then returns to zero when 
the numbers of both species have become equal and reached the value ^N. 

Consider a macroscopic system initially in the pseudo-vacuum of compacts. What happens to the 
entropy if we add equal numbers of particles from two of the four species until capacity is reached? The 
answer is shown in figure[6l At low density, the positional disorder is little affected by the different shapes 
of the particles. The curves overlap close to perfectly. At higher densities the different shapes dictate the 
presence or absence of ordering tendencies and the types of ordering realized. Particles 1,1 or 2,2 or 3,4 
interlink directly. Two close-packed particles 1,2 or 3,3 or 4,4 have two vacuum elements in a particular 
sequence between them. 

With increasing (averaged) densities Ni - N2 or N3 - N4 the entropy varies along the same curve 
even though the associated equilibrium states are very different. In the case of particles 1,2 the shapes 
favor single-species clustering when the space becomes crowded. Mixed-species clustering is favored in 
the case of particles 3,4. When capacity is reached, the entropy has returned to zero. The equilibrium 
state is then fuUy phase separated in one case, consisting of two equal-size single-species clusters, or 
homogeneous in the other case, consisting of one mixed-species cluster. 




Figure 6. Entropy Slk^ per lattice site as function of the density Nm = N^' of two species of compacts 
from table[5] 
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Particles 1,4 (or 2,3) are close-packed with three vacuum elements between them. By contrast, close- 
packed particles 1,1 interlink directly and close-packed particles 4,4 have two vacuum elements between 
them. With increasing (averaged) densities iVi - N4, single-species clustering crowds out mixed-species 
clustering. Capacity is reached earlier. The entropy returns to zero as before. The equilibrium is fully 
phase-separated. The two clusters are of unequal size. 

Particles 1,3 (or 2,4) are separated by only one vacuum element when close-packed. The consequence 
is that with increasing (averaged) densities A^i - N3 single-species and mixed-species clustering are 
equally favorable (compared to loose particles). The system reaches capacity in an amorphous state. The 
entropy stays nonzero. The higher rise of the entropy compared to the previous case is explained by the 
smaller size of close-packed mixed-species pairs, which increases the positional disorder at equilibrium. 



4. Statistical mechanics 

The statistical mechanical analysis of the particles with shapes identified in section [2] can be per- 
formed for open or closed systems. Here we consider an open system. Wu's analysis IIS.1 for a generic 
situation starts from the expression 

Z= E Wi{N,„})eJ-Z'-f^) (5) 



{N,„} 



kBT 



for the grandcanonical partition function, where are the particle energies and W[{Nm]] the multiplic- 
ity function^ That analysis produces the general result 

Z^Yll ^ , (6) 



where the (real, positive) w,„ are the solutions of the coupled nonlinear algebraic equations, 

— = ln(l-H S'm'm In 1^^;— I . (7) 



m' 



The capacity constants Am, and the mutual interaction coefficients gmm' are tabulated in section[2] Arbi- 
trary energies Cm can be assigned to each particle species. With the Wm from 0. the average numbers of 
particles can be derived from (5) via 

51nZ 

{Nm)^-kBT- , (8) 

which, when carried out using (6) and Q, leads to the linear coupled equation^ 



w, 



z{Nm) + Y.gmm'(Nm')^Am. (9) 



In similar fashion we can derive from (6) and l[7) correlations between particle populations, specifically 
the covariances ^ 

{{NmNm')) = {NmNm') - {NmXNm') = [k^T]^ ^ . (10) 

dcmdem' 

The entropy inferred from (6) can be expressed as a function of the {Nm) alone, namely by the function 
S{{{Nm)]) from g). 

It is instructive to compare this method with the transfer matrix method (ll] in the context of the 
Ising chain l[l). The latter operates with coupled degrees of freedom of minimal structure whereas the 
former operates with degrees of freedom that are no longer coupled but have more complex structures. 
Ising spins are tied to lattice sites whereas particles are floating. The number N of Ising spins is fixed 
whereas the numbers Nm of particles from each species are fluctuating. Ising spins at different sites are 



■^The chemical potential is zero for aU species, a consequence of the energy scales chosen in section|2] 

*The analysis assumes that W » 1 and {Nm) » 1. Any part of a nonzero A,,, that is of 0(1) is irrelevant. Wu's derivation of (6) 
proceeds via (9). The calculation of (5) from (8) via (6) and is merely a check of internal consistency. 
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distinguishable whereas particles from the same species are not. The canonical partition function Zm 
from the transfer matrix analysis is related to the grand partition function (6) via 

^^gfipv/fcar^^^ (11) 

where Ep^ is the energy of the pseudo-vacuum in use. 

For a brief demonstration of how particles with interlinking motifs and different energies compete 
for space in the presence of thermal fluctuations and produce long-range ordering tendencies at low T 
as a consequence, we assign the energies of Q) at = as tabulated. All thermodynamic quantities of 
interest can then conveniently be expressed in terms of the functiorjf] 



ulKj, Kl) = coshKf + \l sinh^ Kj + e^J^i , Kj±^—, Kl^^^. (12) 

The famUiar transfer matrix result then reads 

ZN^[u{Kj,KL)e-'^^f . (13) 

4.1. Analysis from Ititi •••>2 

The physically relevant solution of equations l[7) for the four species of hosts and tags introduced in 
section[2A]is 



= If 2 = u{Kj,KL)e^' - 1, w-i- Wi- w\e (14) 
The grand partition function © becomes 

Z^[u{K,,KL)e-^'f , (15) 



consistent with (13) via dill . Equations (9) yield 



<M> = {N2) = <Af3> — , (m ^(Na)^\ . (16) 

In figure [7]we show contour plots of the population densities of tags m - 1,2 (111.111) and hosts 
m - 3,4 (UU. I III)- Lowering T at fixed J,L means moving from the center radially outward in a 
given direction. Any particular population density either increases or decreases from the common value 
{N„j) = I at r = 00. In region $2 all (Nm) decrease and reach zero at T = 0. All particle energies Cm are 
positive. Throughout region <I'4 the host energies are negative. The host population densities (N^) - (N4) 
increase as T is lowered and reach the value ^ at T = 0. In part of this region, at / < 0, the tags have 
negative energies as well, but their energy density is less negative than that of hosts. Hence tags are 
crowded out by hosts. In region <l>i the host population diminishes and the tag population proliferates 
toward (ATi) - (N2) = | as r is lowered. In part of this region, for L > 0, the hosts have lower (negative) 
excitation energies than the tags. Nevertheless, the hosts are crowded out by the tags due to their larger 
size. 

On the (t>i-92 boundary (/ = 0, L < 0) the hosts have em > and the tags Cm = 0. As T is lowered the 
host population decreases, leaving more room for tags and elements of pseudo-vacuum, whose motifs 
do not interlink directly. At low T the few remaining hosts act as surfactants between segments of tags 
inside and segments of pseudo-vacuum outside d]. The ground state is fourfold degenerate. Two vectors 
contain no particles. The other two are solids formed by tags of one or the other kind. The density of tags 
averaged over the four vectors is (Ni) - {N2) - \ ■ 

On the (^2 - boundary (L - J 12 > 0) the particles with Cm > are tags and the particles with e,„ - 
are hosts - a switch with drastic consequences. As T is lowered the tag population decreases, surren- 
dering the lattice to hosts and elements of pseudo-vacuum, whose motifs do interlink directly. A highly 



^The mapping which relates the special cases h = and I = of Hamiltonian (T) with open boundary conditions by simple 
parameter transcription J -h, J 120] is evident in the partition function but has no impact on what is demonstrated here. 
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Figure 7. (Color online) Average number (per site) (a) (Ni) = <Af2> of tags and (b) (N^) = (N4) of hosts 
versus Kj,Ki. Plot (a) also represents the population density of FM particles discussed in section |4!2] 
The population density of AFM particles from that same set is the left-right mirror image of plot (a). The 
dashed lines indicate the r = phase boundaries of figure[l](b). 

degenerate ground state ensues. The limiting host population densities (16) are ( - {N4) - (5+ \/5)~^ - 
0.138. 

Finally, on the - (^4 boundary (L = - J/2 > 0), all particles have < 0. Hosts and tags have equal 
negative energy densities. All configurations of hosts and tags covering the lattice solidly have the same 
energy. The resultant entropy is the same as that along the <^2 - ^4 boundary. The host population density 
remains the same, - (N4) - [5+ v^)~^ - 0.138. The tag population density, calculated from (16), is 
(Ni) = im = (2 ^ 0.224. 



4.2. Analysis from I nUII •••>4 

The grand partition function for the four compacts introduced in section lOl becomes 

again consistent with (13) via (TT) . and where 

wi- wz- u{Kj,KL)e^' - 1, w-i-W4-[\ + wi)e~^^' - 1, 
are the solutions of equations (7), and equations (9) yield 

W-i - - \ W\ 



1 



{Ni) = {N2) = - 



2 2w\W-i+ W\+ w-i 



2 2 Wi If 3 + If 1 + If 3 



(17) 



(18) 



(19) 



The population densities of the FM particles 1,2 (| fl, 1 11) and the AFM particles 3,4 (| If, | U) have a 
mirror-image relationship (figure[7)H Lowering T from infinity either leads to a depletion of FM particles 
(region O2) or to a depletion of AFM particles (region ti) or to a depletion of both kinds of particles 
(region <1>4). 

If T is lowered on the - ^2 border, all four particles have equal and negative energies. What mat- 
ters for the nature of the ground state is that FM particles and AFM particles interlink among themselves 
(albeit in different manner) but not with each other. The ground state thus consists of the four states that 
are either packed with FM particles or with AFM particles. 



^Tag 2 from tableland FM compact 2 from table|5]have identical motifs and, therefore, identical (absolute) energies. However, 
they are different particles existing in different worlds. Nevertheless, their population densities are identical ath = 0. 
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Along the Oi -<I>4 border, FM particles have zero energy and AFM particles have positive energies. The 
lowest energy level is highly degenerate, comprising states populated by FM particles in all permissible 
configurations. The average population density of these compacts is (Ni) - (N2) - (2 v^)"^ - 0.224. The 
role of the AFM particles along the ^2 - ^4 border is the same as that of the FM particles along the - O4 
border in all respects except the way the motifs interlink. 

4.3. Analysis from | If • ••>i 

Equations of the five species of hosts, tags, and hybrids introduced in section l2^ have the solutions 



W3^u{Kj,KL)e'^' -1, 1+ ^^-4Kj^ 1+ W5 ^^-4Kl^ (20a) 

(1 + w^y W3W4 

wi = e"^^S W2 = — , (20b) 

1+ W4 Wi If 3 - 1 

implying 

Z^[u{Kj,KL)e'^'f . (21) 

The solutions of equations (9) are again straightforward but, in this instance, unwieldy. Symmetries 
dictate that (7V2> - (Ns). The population densities as presented graphically in figure|8]are more complex. 
Two of the species play only minor parts in the statistical mechanics. Hosts 2 (If i ff) and hybrids 5 (i i f ii) 
are not abundantly present for any combination of coupling constants or any temperature. They barely 
make it to an average of five particles per hundred lattice sites under the most favorable equilibrium 
circumstances, whereas hosts 1 (If I U I) make it to an average of up to 25 and the tags 3 (i U), 4 (fi fi 
, if it) to averages of up to 50. 

At L > 1/1/2, host 1 has the lowest energy of all by a considerable margin. Its population density is 
dominant in that region. At / > 0, L < the equilibrium state is dominated by tags 4 and at / < 0, L < by 
tags 3. In both regions these low-energy tags can only exist as attachments to higher-energy hosts. As T is 
lowered fewer and fewer hosts have more and more tags attached to them. At T = the last host particle 
is replaced by a low-energy tag in a configuration where each tag is attached to its neighbor on the left 
around the chain. 

In the remaining two regions at < L < |/|/2 two or three species of particles are in competition 
for dominance with energy, size, and category all factoring in. The smallest size among the low-energy 
particles happens to win in both regions as the temperature is turned down: tag 3 at / < and tag 4 at 
/ > 0. There is no region where the total particle population is thinned out entirely as T ^ 0. Tags 3 have 
zero energy. They only disappear if they are crowded out by particles with negative energy. Our choice 
of a pseudo-vacuum that does not coincide with the physical vacuum. 

Along the <l>i - <I>2 border tags 3,4 both have zero energy while aU other particles have positive en- 
ergies. Tags 3 do not interlink with tags 4. The ground-state degeneracy remains low (fourfold). On the 
$1 - (I>4 border we have again two zero-energy species (hosts 1 and tags 3) and three positive-energy 
species. However, the former two do interlink. The ground-state degeneracy is high. The population den- 
sities are (Ni) = (5-h \/5)~'^ 0.138 for the hosts and {N3) = (2 v^)"^ ^ 0.224 for the tags. 

A considerably more complex scenario unfolds on the <1>2-'1'4 border. Four species of particles of three 
different sizes and with three different negative energies are in competition. Only tag 3 (with zero energy) 
stands by. Of the four competing particles, 1,2,5 can host and 4,5 can be hosted. Hosts 2 and hybrids 5 are 
absent at T = except on the phase boundary. There we have (Afi> - {N4) - | , {N2) - (Nb) = ^ . At T > 
the particles 2,5 maintain significant populations only near the phase boundary 



4.4. Entropy 

The three sets of particles produce distinct entropy landscapes in the space of population densities as 
discussed in sectionjS] This diversity is attributable to the different structures of the particle species and 
their mutual statistical interactions as encapsulated in the entropy expression (4). However, when applied 
to a particular physical context, where aU particle energies are functions of Hamiltonian parameters. 
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the function S{{{Nm)\) inferred from jiTI encodes, at any given temperature, the same entropy landscape 
for all three sets of particles in the space of these parameters. That landscape is shown in figure|8](d) for 
the situation at hand. 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



Kj Kj 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



Kj Kj 

Figure 8. (Color online) Average number (per site) (a) (iVi), (b) {N2) = (Ns), and (c) (^4) versus scaled 
variables Kj,Ki. The corresponding plot for {N3) is as shown in figure [7] (a). The graphs of {N3) and 
{Ni) are not mirror images of each other. The graph of (iVi) + {N2} is identical to the graph shown in 
figure[7](b). Panel (d) shows the entropy Slk^ per lattice site as function of the scaled variables Kj,Ki . 



For given / and L the entropy is a smooth and monotonously increasing function of T. Magnetic short- 
range order (one of three kinds) establishes itself gradually and turns into long-range order at T = 0. The 
thermal fluctuations are stronger at L > than at L < if / owing to competing nn and nnn couplings. 
S/fcfi decreasing smoothly from ln2 0.693 at T = cxd to zero at T = for all parameter combinations 
except I/I - 2L. Here S/ks approaches ln([l -1- ^5] /2) ^ 0.481 in the low-T limit. The implied ground-state 
degeneracy has a natural yet different interpretation in the context of each set of statistically interacting 
particles as shown. 



^Alternatively derived from (13) via S = tedrlnZiv/dr. 
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5. Structures from interactions 



In this work we have investigated structures caused by interactions in different ways. On a small 
scale, particles of various shapes are assembled from building blocks with minimal structure (Ising spins) 
by the nn and nnn couplings of Hamiltonian (1). These particles with structures over the range of a few 
lattice sites are free of any inter-particle binding forces, not merely at low density but also at high density. 
Nevertheless, by virtue of their specific shapes, these particles are apt to assemble further structures on 
intermediate and large scales in the form of positional ordering of various kinds. 

On an intermediate scale we have described the assembly of clusters of compacts, a process driven 
entirely by shapes and limited space. We have also described host particles consisting of two amphiphilic 
parts bracketing uniform arrays of tags or disordered mixtures of tags and hybrids, thus assuming the 
role of a surfactant or a membrane. Host particles are akin to micelles in this context. On a larger scale we 
have described the formation of macroscopically ordered patterns of particles from one or two species 
in a crowded environment. The type of ordering is determined solely by the energies, shapes, and sizes 
of the particles. 

The structures stabilized on the smallest scale by the coupled Ising spins depend on our choice of 
reference state (pseudo-vacuum), hence the three sets of particle species in tables [l] [3] [5] However, in 
any given equilibrium state, the same large-scale ordering tendencies are produced by particles from 
different sets, i.e. by particles with different shapes, sizes, and energies. 

The methodology developed here and in reference |8] has natural applications in research areas of 
strong current interest including jamming of granular matter in narrow channels jl^l and DNA over- 
stretching I.l8i.il9i1 . 
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B3aCM03B'}13aHi MOTMBM 1 yiaHALUa())TI/l eHTpOnil CTaTI/ICTI/IHHO 
B3aCM0AiK)m/IX MaCTI/IHOK 

n. A- -Hi^, r. Mto/i/iepP, M. Kap6a>P 

cPisnHHuPi (|)aKy;ibTeT, YHiBepcuTeT Pofl-AR/ieHfly, KiHrcTOH Pofl-Afi/ieHfl 02881, CLLIA 
tpisnHHi/iti c|)aKy;ibTeT, ripHUHnCi yHiBepcuTer Bynnepra/ifl, 42097 Bynneprayib, HiMeH4HHa 

BHK0pHCT0ByK)4H 5 = 1/2 ISHHriBCbKHi^ ;iaHL4K)>K0K 3 OflHOpiflHOK) BSaCMOflieK) Hati6;iM>KHHX i HacTynHMX 3a 

HaPi6;inxHHMW cycifliB, no6yflOBaHO cwcreMy HesacfiiKCOBaHnx nacTWHOK, akI xapaKiepusyKDibCfl MOTWBaMU lue- 
CTH noc/iiflOBHMX /lOKa/ibHux cnihiB. BsacMOflifl cniHlB cnpi/i4HHfle rpynyBaHHfl nacTHHOK, akI, b cbokd nepry, 
He sa/iexaib Bifl eneprii^ BsacMOflii HaBiib npn bucokIii rycTWHi. Bci MiKpocraHU yTBOpeni KOHcfiirypai^iaMW 

HaCTWHOK 3 OflHOrO i3 TpbOX pi3HWX Ha6opiB, flKi BiflnOBiflaKDTb 36yfl>KeHH51M nceBflOBaKyyMiB, LJ40 n0B'a3aHi 

i3 OCHOBHWMH CTaHaMW nepioflUHHOCTi OflHH, p,Ba i HOTupw. Mothbh 4acTHH0K i eyieweHTH nceBflOBaKyyMy 
o6'eflHyKiTbca b flBOx cniyibHi/ix By3;i0Bnx 3MiHHHx. CraTHCTUHHa B3aeM0flia Mix nacTUHKaMi/i e 3aK0A0BaH0K) 
B y3ara;ibHeH0My npuHLii/ini flay/ii, mo oni/icye m po3MimeHHfl oflHOi nacTUHKn SMlHioe MO>Kyii/iBOCTi a/ia po3- 
MimeHHfl nofla/ibiuMx nacTi/iHOK. B CTari/iCTUHHO MexaHiHHOMy aHa;ii3i flOBiyibHi eHeprii Moxyib craBi/iTucb y 

BiflnOBiflHiCTb BCiM COpiaM HaCTWHOK. EHTpOnifl £ (fiyHKL^iCK) 3aCe;ieH0CTi HaCTWHOK. OC06/1HBOCTi CTaTHCTWHHOI 

B3aeM0fliT npo3opo B6yflOBaHi b u,e\/\ Bwpa3. EHeprii i crpyKTypn nacTHHOK bhk;ik)4ho BH3Ha4aK)Tb Bnop^flKyBa- 
Hha npn HM3bKnx TeMnepaiypax. 3a oco6;ihbwx yMOB HacTWHKW Moxyib 6yTM 3aMiHeHi cfiyHflaMeHTa/ibHiiunMW 
4acTHHKaMM 3 KopoTLUHMW MOTMBaMM, w,o B3aeM03B'fl3yKiTbca OflHicK) cniyibHOK) By3;iOBOio smIhhok). 3yM0B- 
;ieHi B3aeM0flieK) crpyKiypw BWHWKaiOTb Ha flBOx piBH^x: nacTUHKH 3 c|)opMaMM, akI yiBopehi B3aeM0fliK)4HMH 
cninaMU, i TeHfleHL4ii' flo fla/ieKOcaxHoro BnopflflKyBahHfl 4acTHH0K 3 ((lopMaMU, akI CTarncTUHHO B3aeM0fliK)Tb. 

KniOMOBi c;iOBa: npi/iHu,i/in flayjii, ^acTmnm 3 ipopMaMi/i, flpodosa CTaTi/icmKa, Mo^enb l3MHra, conhoHi/i 
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